Small x nonlinear evolution with impact parameter and the structure function data 
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The nonlinear Balitsky-Kovchegov equation at small x is solved numerically, incorporating impact 
parameter dependence. Confinement is modeled by including effective gluon mass in the dipole 
evolution kernel, which regulates the splitting of dipoles with large sizes. It is shown, that the 
solution is sensitive to different implementations of the mass in the kernel. In addition, running 
coupling effects are taken into account in this analysis. Finally, a comparison of the calculations 
using the dipole framework with the inclusive data from HERA on the structure functions F2 and 
Fl is performed. 



I. INTRODUCTION 

With particle colliders extending the energy frontier, the need to understand QCD in the high energy limit becomes 
essential. The new kinematic regime explored at the LHC and potentially in the future deep inelastic scattering (DIS) 
machines, EIC [T], Q and LHeC p| H[, requires detailed analyses of the high energy and density limit in QCD. At 
these high energies it is expected that the parton densities will become very large. In particular, as the Bjorken x 
becomes very small one needs to take into account the effects of parton saturation p| [f| . Evolution into the small 
x region, which incorporates parton saturation phenomenon, is governed by the nonlinear Balitsky-Kovchegov (BK) 
[7H12j equation. A general framework which systematically incorporates the effects of the large gluon density is the 
Color Glass Condensate (CGC) [l3|, [l4[ model. Within the CGC model the evolution into the small x region is given 
by the renormalization group - type equation, the JIMWLK equation [loH^H • The CGC framework contains the BK 
evolution equation as well as the evolution of the higher point gluon correlators, for a recent study see [22| . The BK 
equation is an extension of the linear BFKL evolution equation for small x [23l425| as it takes into account parton 
recombination effects. These effects are included through the additional nonlinear term in parton density. As a result 
the solution to this equation cannot exceed unity, which is the unitarity bound for the dipole (i.e. quark-antiquark 
pair)-target scattering amplitude. 

Despite the fact that the BK equation is closed and a relatively simple nondinear equation, no exact analytical 
solution yet exists. Nevertheless, there have been numerous analytical studies [26l - |29| as well as extensive numerical 
analyses [30M36| , and the properties of the solution are currently very well known. The solution has also been used to 
successfully describe the experimental data on the structure function F 2 [3l|, H3, HH . Furthermore, it has been used 
in the prediction of a large number of processes in hadron and heavy ion collisions, such as multiplicites and single 
inclusive spectra [39l44ll ]. 

In most of these analyses one utilizes an assumption of the impact parameter independence of the solution. To be 
precise, the dipole target scattering amplitude in this approximation depends only on the dipole size and not on the 
position in impact parameter space. This leads to a significant simplification of the problem and drastically reduces 
the CPU time needed to numerically solve the equation. On the other hand, the impact parameter is an important 
ingredient for many of the phenomenological predictions. For example, the total multiplicities in heavy ion collisions 
depend strongly on the ccntrality and hence knowledge of the impact parameter distribution of the partons is essential. 
Usually this problem is circumvented by assuming an average parton density distribution for each of the centralities, 
(for a more refined approach which includes nucleon configurations in the nucleus in the initial condition see [42j). 
In the context of the saturation physics it is expected that the saturation scale Q s , which characterizes the dense 
system, is impact parameter dependent (i.e. Q s (b)). The saturation scale will thus have a larger value close to the 
dense center of the interaction region and smaller value at the periphery of the interaction region where the partons 
form a dilute system. The knowledge of the parton density distribution in impact parameter is essential not only 
in the context of heavy ion collisions but also for hadronic collisions (for example in the problem of the underlying 
event) and in particular for diffractive processes. For example, in exclusive diffractive production of vector mesons in 
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DIS the momentum transfer dependence crucially depends on the impact parameter profile of the dipole scattering 
amplitude 

The impact parameter dependence has been taken into account in Monte Carlo simulations based on dipole evolution 
and scattering [461 - 151} as well as in numerical solutions to the BK evolution equation, see [lH, [53| and [54| . It has been 
also discussed in the context of the conformal properties of the equation [55j. The results so far indicate that there 
are important modifications to the solution when the impact parameter is taken into consideration. In particular, it 
has been found that the parton (or more precisely dipole) density distribution in impact parameter possesses long 

range Coulomb-like power tails m m m m, which are the direct consequence of the form of the perturbative 

branching kernel. These tails need to be regularized by the appropriate cuts on the large dipole sizes which mimic 
confinement effects [4H - [5ll IBlj . The other important observation was that the dipole scattering amplitude decreases 
with increasing dipole size at large dipole sizes (when it is evaluated at fixed value of impact parameter). This has 
to be contrasted with the impact parameter independent solutions for which the amplitude always saturates to unity 
for arbitrarily large values of the dipole size. 

In the earlier work (53j we explored the dynamics of the BK equation with impact parameter by taking into account 
subleading effects such as kincmatical cuts and the running of the strong coupling. The goal of this paper is to use 
the solutions to the BK equation with impact parameter dependence to compute the cross sections and structure 
functions for the deep inelastic process. In order to perform this analysis we introduce the initial condition with 
physical scales and also modify the branching kernel to account for non-perturbative confinement effects. This is done 
by introducing a mass parameter m into the kernel which modifies the long distance behavior of the dipole-target 
scattering amplitude. This parameter restricts the splitting of dipoles into daughter dipoles which are larger than 
^max = — in order to account for confinement. With such a setup we then compute Fn and Fl structure functions 
using the resulting solutions and compare them with experimental data from HERA [56l . . 

The resulting dipole scattering amplitude was then compared with the parametrizations available in the literature 
PH which include impact parameter. In particular, we find that although the dynamically generated amplitude from 
the BK equation is similar for small values of the dipole size to the Glauber-Mueller like parametrization, for larger 
dipole sizes one observes notable differences. To be precise, we find that the BK equation generates solutions which 
possess specific correlations between the dipole size and impact parameter, an effect which is totally absent in the 
Glauber-Mueller type parametrizations. 

The outline of the paper is the following: in Sec. [IT] we state the basic formulae for the inclusive cross section within 
the framework of the dipole model. In Sec. IHII we briefly discuss the most substantial features of the solution with the 
impact parameter dependence and discuss the differences with respect to the impact parameter independent scenario. 
We then introduce modifications to the evolution kernel by including the mass which mimics confinement. We discuss 
the properties of the solutions which result from these modifications. Both fixed coupling and running coupling cases 
are considered as well as the various methods that we used to implement the mass parameter which regulates the 
large dipole sizes. Later, in Sec. II VI we show the results for the dipole amplitude and we make the first comparison to 
the data for F 2 and Fl. Finally, in Sec|V]we state the conclusions. 



II. DIS INCLUSIVE CROSS SECTION WITHIN THE DIPOLE MODEL 

The dipole model [H, [59| is a very useful tool in evaluating many processes at small values of x. One of the 
advantages of this approach is the possibility of including multiple parton scattering effects. It has been originally 
formulated for the description of deep inelastic lepton-proton (or nucleus) scattering at small x. In this picture, 
utilizing the leading logarithmic approximation in x, the incoming electron emits a virtual photon which fluctuates 
into a quark-antiquark pair, a dipole. The color dipole then subsequently interacts with the parton constituents of 



the nucleon, as is illustrated in Fig. 1(a) The interaction of the dipole pair with the target is given by the scattering 
amplitude N. The qq pair is characterized by a dipole size which is defined as a separation distance of the color charges 
xoi = xo — xi (where Xo and Xi are the positions of the q and q in transverse space). 1 The transverse momentum of 
the quarks in the dipole is of the order of ~ where large dipoles correspond to the infra-red region and need to be 
regulated, as will be discussed in detail in later sections. The interaction of the dipole with the target is described by 
the scattering amplitude N(r, b; Y) which contains all the information about the dynamics of the strong interaction. 
In the following analysis the full dependence of the scattering amplitude on the impact parameter b will be taken 
into account. The evolution of the amplitude in rapidity Y can be represented as emissions of daughter dipoles from 



In this paper we shall denote vector quantities in bold, otherwise they should be read as magnitudes of the associated vector. Also, 
alternatively we will be also using here the notation for the dipole size to be r = xqi and impact parameter b = 
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FIG. 1: Plot (a): schematic representation of the dipole picture in DIS. The incoming virtual photon 7* splits into a color 
dipole (quark-antiquark pair) of size r which subsequently interacts with a target, where N is the dipole-target scattering 
amplitude. Plot (b): depiction in transverse space of single branching of the parent dipole 01 into two daughter dipoles. The 
sizes of dipoles are denoted by Xij. The impact parameter variables bij of all the dipoles are shown relative to the target. 



the original parent dipole. When the original dipole 01 splits into two dipoles 02 and 12 a new coordinate appears, 
X2. These two daughter dipoles are produced with sizes X12 and X02 at impact parameters t>i2 and t>02 as illustrated 
in Fig. |l(b)| 



The dipole-target amplitude N(r, b; V) at high values of rapidity Y (or small x) is found from the solution to the 
BK evolution equation which can be represented in the following form: 

cWxnx, f d 2 x 2 



'xpx i 

dY 



/d?"x. 
—^K.(x \,xi2,x Q 2\a s ,m) [iV XoX2 + iV X2Xl - iV XnXl - N XoX2 N X2Xl ] . (1) 



In the above equation we used the shorthand notation for the arguments of the amplitude N XiXj = 
N(rij = Xj — Xj, by = i(x, +Xj); Y) which depends on the two transverse positions Xj and Xj and on the rapidity Y. 
The branching kernel )C(xoi , £12, ^02; ce s , m) depends on the dipole sizes involved and contains all information about 
the splitting of the dipoles. In addition it depends on the running coupling a s . We have also indicated that it depends 
on the infra-red cutoff m imposed on large dipoles. 

The solution to Eq. ([T]) is the dipole-target scattering amplitude for arbitrarily small x. In order to compute the 
structure functions F 2 and Fl for the proton we use the following standard formulae in the dipole picture in the 
transverse coordinate representation 

F 2 (Q 2 , x) = — ^— f d 2 r [ dz (|* T (r, z, Q 2 )\ 2 + |* £ (r, z, Q 2 )\ 2 ) a dip (r, x) , (2) 

and 

F L (Q 2 ,x) = - ^ [ d 2 r [ dz\^f L (r,z,Q 2 )\ 2 (T dip (r,x) . (3) 

Here, adip is the (dimensionful) dipole cross section obtained from the (dimensionless) scattering amplitude by inte- 
grating over the impact parameter 

a dip (r, x) = 2 J d 2 b N(r, b; Y) , Y = In 1/x . (4) 

The ^(r, Q 2 , Y)t/l functions are the photon wave functions. They describe the dissociation of a photon into a qq 
pair and can be calculated from perturbation theory. The photon wave function has the following form for the case 
of transverse photon polarization 

\* T (r,z,Q 2 )\ 2 = 3 -^J2 e2 f(l z2 + ( 1 - z f] Q) K l (Qfr)+m 2 K 2 (Q f r)) , (5) 
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and for longitudinal polarization 

|^(r^,Q 2 )| 2 = ^E e /( 4 ^ 2 (l-^o 2 (^)) ■ (6) 
/ 

In the above equations Q 2 = z(l — z)Q 2 + mjp, where — Q 2 is the photon virtuality and z, (1 — z) are the fractions 
of the longitudinal momentum of the photon carried by the quarks. In addition -Ko,i are modified Bessel functions of 
the second kind. The summations are over the active quark flavors / of charge e/ and mass to/. 



III. DIPOLE EVOLUTION WITH IMPACT PARAMETER DEPENDENCE 



The solution to the BK equation with impact parameter dependence is found numerically, the description of the 
procedure was outlined in [52l |53| . The technical complication when trying to solve BK with impact parameter is 
the increased number of arguments in the dipole amplitude. In the 6-independent scenario the amplitude depends 
only on two variables: rapidity and dipole size. In the 6-dependent case there are 5 variables: rapidity, dipole size 
(vector, 2-dim.) and impact parameter (vector, also 2-dim.). Alternatively, one can choose the coordinate variables 
to be parametrized by the dipole size, impact parameter and two angles: one parametrizing the absolute orientation 
of the dipole-target system in the coordinate space and the second describing the relative orientation of the dipole 
with respect to the target. With the assumption of the global rotational symmetry (i.e. independence of the first 
angle) the number of variables reduces to 4: rapidity, dipole size, impact parameter and the angle between r and b. 
Such large number of variables requires working with a very large multidimensional grid and it leads to a significant 
increase of computation time per each step of evolution in rapidity. The BK equation was solved numerically by 
discretizing the scattering amplitude in terms of variables (log 10 r, log 10 b, cos 9), where 9 is the angle between impact 
parameter b and dipole size r. The amplitude N(r, b, cos9) was placed on a grid with dimensions 200 r x 200& x 20g. 
More details can be found in Refs. (52l l53j. 

Let us briefly summarize the most important results of (52l . [53j which pertain to the properties of the solution with 
the impact parameter and the differences with respect to the impact parameter independent approximation. We note 
that for the solutions presented in Fig. [5] we used the same initial condition at Y — as in Ref. [Hj], which was taken 
to be 

N° = 1 - exp(- Cl r 2 exp(-c 2 6 2 )) , (7) 

with ci = 10, C2 = 0.5. 

The most distinctive feature of the solution with impact parameter dependence is that the amplitude for large 



dipole sizes goes to zero. This is clearly illustrated in Fig. 2(a) where the solution for fixed value of impact parameter 



b is shown as a function of the dipole size r. This property of the amplitude has to be contrasted with the impact 
parameter independent solution, in which case the amplitude is always equal to unity for sufficiently large dipole sizes. 
The fact that the 6-dependent amplitude drops for large dipole sizes has a rather simple physical interpretation: the 
interaction region possesses finite extension in impact parameter space. The size of this region is set at low rapidity 
by the initial conditions (in our case it is parameter c 2 in Eq. [7]) and is later increased in the course of the evolution 
by the diffusion of the dipoles in transverse coordinate space. The probability of the scattering for dipoles which have 
sizes larger than the extension of the interaction region is very small and therefore the amplitude will tend to go to 
zero for such configurations. The dipole-target amplitude therefore is largest for the scattering of dipoles with sizes 
comparable with the typical size of the target. As a result, the amplitude decreases either for small or for very large 
dipole sizes, which in each case are very different than the extension of the target, yielding a maximum contribution for 
some intermediate sizes of dipoles. The configurations for which the amplitude is small are schematically illustrated 



in Fig. 02 and they correspond to the tails of the distribution depicted in Fig. 2(a) In the course of the evolution in 



rapidity this distribution in the dipole size broadens due to diffusion. As a consequence, the amplitude in the impact 



parameter dependent scenario has two fronts as is evident from Fig. 2(a) The first front (for small dipoles) is similar 
to the front in the impact parameter independent case, and one can define the saturation scale which divides the 
dense (i.e. where N ~ 1 ) and dilute (JVcl) regime for small dipoles. This saturation scale can be parametrized in 
the form Q 2 sL {Y 1 b) = Q\ sL exp(A s iF) with X s l — 4.4 which is consistent with the analytical predictions and with the 
solutions for the impact parameter independent case. The second front expands towards larger dipoles with increasing 
rapidity. One can also define the corresponding 'saturation scale' for large dipoles which can be similarly parametrized 
as Ql R {Y,b) — Qq sR exp(— \ s rY) with \ s r ~ 5.8. Note the '-' sign in the exponent, which originates from the fact 
that the second front expands towards larger dipoles with increasing rapidity, and therefore this second 'saturation 
scale' decreases with rapidity. One has to stress however that the region of large dipoles is going to be heavily modified 



5 




FIG. 2: The dipole scattering amplitude N(r, b; Y) as a function of the dipole size r from the solution to the LL BK equation 
with impact parameter dependence. The strong coupling is fixed a„ = 0.2. The consecutive curves shown in plots are for 
rapidities Y = 10, 20, 30, 40, 50 on the plot (a) and in rapidity intervals of 5, until Y = 50 on plot (b). The dotted-dashed line 
in the left plot is the initial condition at Y — given by Eq. [7] The initial condition is not visible on the right plot as it is very 
close to zero. The orientation of the dipole with respect to the target is such that rib. 



by the non-perturbative effects (see the discussion later in this section) and therefore this second saturation scale is 
most likely only an artefact of the perturbative expansion in the leading logarithmic (LL) approximation. 

This novel feature of the solution in impact parameter dependent case, namely the decrease of the amplitude at 
large dipole sizes is directly related to the profile in the impact parameter space. It was found in the LL case l53j 
that the solution has a power-like tail in impact parameter. This is related to the fact that there are no mass scales 
in the perturbative evolution and hence the interaction is long range [60l IfjiT ] . As we will discuss later in this section, 
the branching kernel in the equation needs to be modified by including the effective gluon mass in order to regulate 
the power- like behavior of the amplitude for large dipole sizes and include the effects of confinement. 

Another distinctive feature of the solutions with impact parameter dependence is the presence of the strong correla- 
tions between the dipole size and the impact parameter. It has been observed that the amplitude is largest for specific 
configurations and orientations of the dipole size r and impact parameter b vectors. In particular, the amplitude has 



a peak when the dipole size is equal twice the impact parameter. This is clearly illustrated in Fig. 2(b) where the 
peak in the amplitude at r = 2b emerges. In this case there is also a non- negligible dependence on the angle between 
the dipole size vector and the impact parameter. It turns out that the amplitude is largest when the angle between 
the dipole size vector and the impact parameter vector is equal to or tt, that is when the vectors r and b are parallel 
or anti-parallel. In this configuration one of the color charges scatters off the center of the target. This enhancement 
can also be seen analytically from the conformal eigenfunction representation as discussed in [531 ] . 

By integrating the resulting amplitude N(r, b; Y) over the impact parameter b, as in Eq. [H one obtains the dipole 
cross section as a function of the dipole size and rapidity. Even thou gh t he amplitude N never exceeds unity, the 
resulting dipole cross section increases very strongly with the rapidity [60L l6lj. due to the rapid diffusion of dipoles 
in the impact parameter space. In the leading logarithmic approximation this increase is exponential in rapidity, 
fdip ~ exp(Asy), with \b — 2.6 [13, HH]. This behavior is very different from the features observed in the b- 
independent solution. In the latter case one assumes that the dipole amplitude does not depend on the position in 
coordinate space, but only on the absolute value of the dipole size and rapidity, N(r; Y). The solution still saturates 
to unity which is the fixed point of the equation in this approximation as well. To obtain the dipole cross section one 
needs then to multiply the amplitude by a dimensionful coefficient i.e. 

ffdip = So N(r; Y) , 

where So can be interpreted as the integral 

S = [ d 2 b 
Jtz 

over the interaction region 71 in the impact parameter and is entirely introduced by hand. The behavior of So on 
rapidity (whether it is constant or increasing) is thus not determined by the 6-independent BK evolution equation. 
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(a) Small dipole (b) Large dipole scattering, 

scattering. 



FIG. 3: Two different dipole-target configurations for which the dipole scattering amplitude is small (away from unitarity 
limit). 



We see therefore, that the 6-dependent solution to the BK equation, when integrated over the impact parameter 
does not reduce to the solution in the previously studied approximation when the impact parameter is neglected. This 
is an important qualitative and quantitative difference between these two solutions. In fact, the solution with impact 
parameter dependence is more physically motivated as it gives the increase of the dipole cross section with rapidity. 
It will also naturally lead to the increase with the energy of the diffractive slope for the vector meson production. We 
note however, that the results in the LL approximation, even with the running coupling, are not compatible with the 
experimental data. This is due to the fact that the diffusion in impact parameter is very strong (since it is driven 
by the purely perturbative physics) and results in a fast increase of the interaction radius which is not supported by 
the data. In order to tame this growth further corrections are necessary, in particular the inclusion of subleading 
corrections and a mass parameter into the evolution as will be discussed later in this section. 



A. Initial condition for the evolution towards small x 



For the rest of the numerical simulations presented in this paper we will use the initial conditions of the similar form 
as in ([7]) but with parameters which were adjusted to obtain the predictions consistent with experimental data. We 
will use as our initial condition the parametrization from Ref. [44| , where the parametrization in the Glauber-Mueller 
form was used 



N GM (r, b; Y = In l/x) = 1 - cxp \-—r*xg(x, r] 2 )T(b) j , (8) 

with 

T(b) = ^e*€ . (9) 

The formula is used as an initial condition for the BK evolution for dipoles with sizes smaller than the cutoff 
1/m. For dipoles larger than the cutoff the initial condition is set to zero ( see the discussion in the next subsection). 
In formula ((5J) the function xg(x,r) 2 ) is the integrated gluon density function and T(b) is the density profile of the 
target in transverse space with Bq = 4 GeV -2 . This parameter was set to fit the t slope of the diffractive J/^S 
production in (45|. Also, the scale rj was set according to [45| to be equal to rj 1 = Q + ?y 2 with parameters C = 4 and 
t]q = 1.16 GeV 2 . The integrated gluon density in © was also taken from fits performed 45 1. We use © as the initial 



condition at Yq = lnl/a;o, xq = 10 -2 and evolve the amplitude with the BK equation to obtain the solution at lower 
values of x < xq. We also note that the initial condition © depends only on the absolute values of the dipole size 
and impact parameter. The nontrivial dependence on the angle between vectors r and b is not present in the initial 
condition, instead being dynamically generated when the initial condition is evolved with the BK equation. 



B. Including the effective gluon mass into the evolution kernel 

Currently it is unknown how to introduce a massive cutoff on a fundamental level into the small x evolution as 
it is an entirely non-perturbative problem. We have tested various prescriptions and found that there is a rather 
large sensitivity of the resulting solutions to the details of confinement implementation. In addition, there is a strong 
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dependence of the solutions on the way the running coupling is regularized. This sensitivity stems from the behavior 
of the 6-dependent solution at large dipole sizes, as discussed previously. One important difference to note between 
the solution with and without the impact parameter dependence is that in the latter case the running of the strong 
coupling is naturally regularized by the saturation scale, provided the latter is in the semi-hard regime. For example, 
it was observed in Ref. [34j that different prescriptions of the regularizations for the running coupling gave similar 
results in the case of the non-linear evolution. The emergence of the semi-hard saturation scale Q s and its role as an 
infrared cutoff is one of the most prominent and useful features of the nonlinear evolution. In the case when the impact 
parameter is taken into account, the saturation scale strongly varies with b. In the dense region this scale is large, 
and is providing a natural cutoff for the running coupling in the same way as in the impact parameter independent 
solutions. There is however always a peripheral interaction region in impact parameter where the scattering amplitude 
is small and the system is dilute. Consequently the coupling is not regularized in this region by the saturation scale 
which is very small for large b. As a result, the solution in the dilute peripheral regime is governed by the linear 
evolution and becomes extremely sensitive to the region of large dipole sizes. 

The specific value of the massive cutoff which is implemented into the evolution kernel should correspond to the 
non-perturbative scale which is related to confinement. Lattice simulations [62ll63l| suggest the presence of an effective 
gluon mass which would regulate the large-dipole size regime. In the case of BK equation an important ingredient that 
one has to take into account is the fact that the evolution in impact parameter is strongly correlated with the evolution 
of the dipole sizes. Therefore the cutoff on the latter will crucially influence the size of the interaction region in impact 
parameter and its variation with the collision energy. In the Monte Carlo analysis of dipole evolution in l5Ct 16411 the 
parameter r max was set to be around 3 GeV -1 (to be precise it was set to 2.9 GeV -1 in [64| and 3.1 GeV -1 in [50|). We 
set the value of the cutoff to be of the same order, i.e.m = = 0.35 GeV which corresponds to r max ~ 2.86 GeV" 1 . 

Let us finally note here that the impact parameter profile can be accessed through the measurement of the diffractive 
production of the vector mesons [3 EH- From the experimental data 65] at \t\ < 1.2 GeV 2 it is known that the 
diffractive slope B D of the J/* production ~ e Bot ) is of the order of ~ 4.57 GeV~ 2 for Q 2 < 1 GeV 2 and 

~ 3.5 GeV~ 2 for Q 2 > 5 GeV 2 in the energy range 40 < W 1P < 160 GeV. It is also slowly growing with the increasing 
energy W of the j*p system. The value of r max we set in the calculation will strongly influence the variation of the 
width of the impact parameter profile with the energy. Therefore in principle r max can be related to the dependence 
of the Bd with the energy, [gfj. 



1. Regularization of the kernel at leading logarithmic accuracy with fixed coupling 

The regularization of the large dipole sizes in the dipole kernel is a purely non-perturbative effect. One possible 
implementation, which is physically motivated, is to introduce the effective gluon mass m into the propagators. This 
mass will set a correlation length r max = 1/m which will limit the propagation of the strong color force. One can 
then derive the dipole kernel by computing the emission of the gluon from the initial qq dipole 67, 68]. The Fourier 
transform of the momentum space expression of the qqg lightcone density into the coordinate space results in the 
expression which contains the modified Bessel functions instead of powers as in the LL expression. This means that 
instead of the long range Coulomb-type interaction present in the perturbative evolution there is now screened Yukawa 
force with finite range given by r max . The modified branching kernel for dipoles with effective gluon mass m has the 
following form [67], Hg| 



K 2 (mxo2) + K 2 (mxi2) — 2Ki(mxo2)Ki(mxi2) — — 

X02 X12 



(10) 



In the limit when the dipole sizes are small compared to the cutoff, Xij <C ^- , the Bessel functions are approximated as 
K\ (mxij ) ~ rax - . ■ I n this limit the kernel (|10p obviously reduces to the well known expression in the LL approximation 

[67T - I691 ] . On the other hand, the production of large dipoles x^ > ~ is exponentially suppressed. This form of the 
modification of the kernel (fTU|) was also used in the later version of the Monte Carlo simulation for the dipole splitting 
and evolution [641 ] . 

By inspecting expression (|10[) it is clear that this kernel does not vanish completely when one of the dipole sizes is 
larger than 1/m but the other is smaller than 1/m. In other words even for very large parent dipole size the above 
kernel permits the splitting, provided one of the daughter dipoles sizes is below the cutoff. It means that despite the 
presence of the cutoff 1/m there is still a diffusion into the region of arbitrarily large dipole sizes. 

As an alternative to the above scenario we have thus used a second prescription where the splitting is suppressed 
whenever any of the daughter dipole sizes is larger than the cutoff. To be precise, we have tried the second ansatz of 
the form 
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4 2 2 = 



a s N c 



1 



1 



(11) 



The kernel is thus set to zero whenever any of the dipoles is larger than 1/m. This form gives more suppression of 
the dipole splitting than the first modification (|10[) and results in an overall slower evolution in rapidity. 



2. Regularization in the presence of the running coupling 



The implementations of the cutoff presented above can be used with a fixed coupling. The running coupling 
correction brings in an additional complication as the form of the kernel changes rather significantly and the coupling 
is no longer a simple multiplicative factor. The running coupling corrections to the BK evolution have been computed 



in two independent calculations [70| and 7l|. The two schemes differ by the non-trivial subtraction term as was 
demonstrated in detail in (36j that we will not consider in this work. For the purpose of the current work we use the 
scheme of (7(1 



Bal 
rcLL 



1 ( «s(>02) 

x 02 \ a s( x i2) 



- 1 



a s (x 



1 



X X2 \ a s( x Q2) 



- 1 



! -01 



(12) 



The other scheme derived |7l| tends to be more time consuming per one evaluation which in the case of the impact 
parameter dependent simulations with large grids greatly extends the time necessary for the evolution. Therefore for 
purely practical reasons we utilized only Balitsky prescription (|12l) . 

Kernel (fT2"|) reduces to the LL kernel with the strong coupling evaluated at the smallest value of the dipole size, for 
configurations of dipole sizes which are strongly ordered. For example, in the case when xqi ~ £02 3> ^12 we have 



K 



Bal 

rcLL 



1 



0^12) 

a s {x 2 m ) 



1 



1 



N c a s {x\ 2 ) 1 



(13) 



and for xq\ <C £02 ~ X12 we can take 



"3(^2) 



1 and obtain 



K 



Bal 
rcLL 



a s {xl 1 )N c 



Thus for the reference we have also used an alternative prescription which is the minimal dipole scenario defined as 



K?cLL - a s (imn{x 



2 2 

01' X 12i : 



Nr. 



x 01 



TV X 12 X 02 



(14) 



As we shall see later, even though formally kernel (fT2")) reduces to (TT4")) . at least in the cases considered above, there 
are notable numerical differences between the two prescriptions (|12I14[) . In particular, we found that the evolution 
with kernel (TT21 is significantly slower than with flT4l which is consistent with previous numerical results (3(| obtained 
without impact parameter. 

In the case of the kernel with running coupling (|12l) the implementation of the cutoff in the form analogous to /C^£ 
is not entirely trivial due to the rather complicated form of (| 12|) . We have thus used the simplest scenario, imposing 
the cuts on the daughter dipoles as in (fTTj) . To be precise, for the case of the running coupling with scenario [7(| we 
used 



t-Bal 
rcLL,m 



a s {xl l )N c 



a s (x< 



02) 



x 02 \ a s( x 12) 



as( x l 2 ) 

X 12 \ a s( x 02) 



"12-^02 



TTl 



)Q(—-xl 2 ). (15) 



For comparison, we have also used the minimal dipole prescription for the running coupling (|14[) where all of the 
scenarios with masses (llOllip were extended in this case. This is possible as the minimal dipole prescription gives the 
multiplicative factor in the LL kernel, see (fT4")) . 
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FIG. 4: Dipole amplitude as a function of the dipole size from the evolution to the BK equation with mass m = 0.35 GeV. 
The dotted-dashed line is the initial condition |(8j at Y — and each solid line represents a progression in two units of rapidity 
until Y = 8. The solutions are evaluated at fixed impact parameter b — 0.0001 GeV -1 and for the orientation of the dipole 
defined as r _L b. 



In this paper we use the expression of the QCD running coupling with a mass parameter /j, to regulate the strong 
coupling at large dipoles, which is of the following form 2 



a s (x 2 ) 



b In 



A 



QCD 



(16) 



Here b = 12 nf , rif is the number of active flavors. The fj, parameter effectively freezes the coupling at large dipole 

_ 2 m . In our simulations we used a = 0.52 GeV as an infra-red regulator for the strong coupling, 
fi j 



sizes at a s ,fr 
and Aqcd = 



bln[A 

0.25 GeV 



IV. RESULTS 



A. Properties of the dipole amplitude and comparison with the Glauber-Mueller model 

As we saw explicitly in Sec. Mil the solutions to the BK equation with impact parameter dependence possess very 
interesting and novel properties as compared to the case without the impact parameter. In this section we investigate 
in detail the features of the solutions when the mass m is included in the kernel. We study the BK solutions using 
two different prescriptions for the infrared regulators which were introduced in the previous section, see Eqs. (|10[) and 
(fTT|) . We also investigate in detail the differences between the two running coupling scenarios, i.e. (fT4|) and (fT"5|). It 
turns out that the differences between the simulations have rather significant impact on the phenomenology. 

First, we used kernels (I10lll[) . where the running coupling has been implemented using the minimal dipole size as 
the scale (see Eg. 1 14|) . This was done in order to consistently trace the differences between the two implementations 
of the massive cutoff. 



In Fig. 4(a) we present the simulations using kernel (fTU|) . We observe that despite the fact that the effective gluon 
mass m is incorporated into the branching kernel, the scattering amplitude is non- vanishing at arbitrary large dipole 
sizes. For any rapidity there is still a significant diffusion into the large dipole size region. This feature is simple to 
understand by inspecting the form of (|TTI)) . One observes that, for any value of Xoi there are configurations where one 
daughter dipole is very large and above the cutoff but the second daughter dipole can be still below the cutoff j-. 



2 Note that, this is the expression for the running coupling in coordinate space. In the literature one finds various forms of the running 
of the coupling in coordinate space which include different normalizations for the argument, i.e. In -j^-^j with differing values of the 
constant C. We are using here the convention from Ref. \!% where (7 = 1. 
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Impact parameter: 1.000 (GeV" 1 ) | cos(phi): 0.0 | DeltaY: 4.0 | max Y: 8.0 



1.0 - 




Dipole Size (GeV ) 



FIG. 5: Dipole amplitude as a function of the dipole size for the fixed value of impact parameter. Solid line corresponds to the 
simulation with the running coupling using the minimum dipole prescription Eq. (|14p . Dashed line corresponds to the running 
coupling using Eq. (I15|l . In both cases the mass m = 0.35 GeV is included in the kernel, in the form of two theta functions, like 
in Eq. fTTj. Two sets of curves correspond to rapidities Y = 4 and Y = 8. The dashed - dotted curve is the initial condition 
at Y — for both calculations. 



These configurations lead to the non- vanishing contributions of the kernel even at large values of Xo%. Because of this 
effect the evolution with modified kernel (fTTJj) still proceeds into the large dipole regime as indicated by the results in 
Fig.g(g] 



In Fig. |4(b)| we show the solution using the kernel with the theta functions imposed as in (fTTj) . We see from this 
figure that the prescription (fTTj) leads to the solution which completely vanishes for large dipoles in stark contrast 
with the simulation shown in Fig. 4(a) The interesting feature about this scenario is that even though the initial 
condition (shown by dotted-dashed line) was cut at xqi = -jr , the evolution does not respect this cut and moves it to 
a larger value equal to xoi = Jj. This happens because the kernel (fTTj) has cuts only on daughter dipoles and not on 
the parent dipole. It means that the amplitude for sizes of dipoles which are larger than — is non- vanishing. Until 
a point at which the parent dipole is twice the cutoff size there exist configurations where neither emitted daughter 
dipoles are above the cutoff size. These symmetric states exist until xqi = — at which point all of the configurations 
are cut, since at least one daughter dipole is larger than the cutoff. 

From this part of the analysis we conclude that when the massive cutoff is imposed in the form of (fTOj) . the evolution 
still proceeds into the region of large dipole sizes, which is dominated by the large value of the strong coupling constant. 
This has rather important consequences and we found that this scenario is actually disfavoured by the data. 

As a next step, we have performed the comparison of the two solutions using the minimal dipole prescription (|14|) 
and the Balitsky prescription (|15p for the running coupling. The results of this comparison are shown in Fig. [5j In 
both of these cases we have implemented the mass as in (fTTj) . i.e. setting the kernel to zero whenever any of the 
daughter dipoles are larger than the cutoff. It is evident that the kernel with the running coupling given by (TT5l) leads 
to a much slower evolution than the kernel with the running coupling as in (|14j) . For example at rapidity Y = 4, the 
evolution front in the minimal dipole scenario is almost at the same position (for small dipoles) as for the scenario 
(TT5l) at rapidity Y = 8. As we will see later, this will translate into large differences for the observable structure 
functions depending on the prescription used. It has a large effect, in particular on the x slope of F%. This feature 
has been found also in earlier calculation which did not include the impact parameter dependence and masses, see 
Ref. II. 

We have also compared the solutions to the BK equation with the Glauber-Mueller model, Eq. f8j In the latter 
model, the parameters were obtained from a fit to the HERA data [45[ . The initial condition for the solution to BK 
was also taken to be of the form of Eq. ([S]) at Y — which we choose to correspond to x = 0.01 in this calculation. For 
consistency the initial condition for the BK equation is set to zero for dipoles which exceed the cutoff. The comparison 
between the solution to BK and Glauber-Mueller model is presented in Fig. 6(a) We see that the solution to the BK 
equation agrees quite well with the parametrization (|5]) for small values of dipole sizes. On the other hand there are 
sizeable differences in the larger dipole regime, where one sees that (f8j) extends indefinitely whereas the BK solution 
is cut off by the massive regulator. This has a non-negligible impact on the phenomenology of F2 at low Q 2 as will 
be illustrated in the next subsection. 



11 




FIG. 6: Dipole scattering amplitude as a function of the dipole size for the fixed value of impact parameter b. Solid line 
corresponds to the model © with parameters from [4^|. The dashed line is the solution to the BK equation with the kernel 
(fT5l> . The dashed-dotted line denotes the model ^ at xo = 0.01. The cutoff at r max = — (plot (a)) and r max = ~ (plot (b)) 
is also indicated. Impact parameter is fixed to b = 1 GeV -1 . 



The difference in the solutions was examined, when the initial condition is not cut at ^ but at ^ and still the kernel 
(fT"5|) is used with the cutoff ^ in the evolution. The result is shown in Fig. |6(b)| In this case it is evident that the 
cutoff is not moved due to the reasons described earlier in this section. However, there is a peculiar structure of the 
solution, where a second peak of the amplitude emerges for dipoles somewhat smaller than the cutoff. The solution 
in Fig. |6(b)| is nevertheless very close to the solution shown in Fig. 6(a) for small values of dipole sizes, i.e. smaller 
than — . We also note that for the values of the cutoff m used here, which are motivated by the profile in impact 
parameter, both the initial condition and the BK solutions do not completely saturate, although the amplitude is 
close to 1 for large rapidities. 

It should be also noted that the solutions to the BK equation with the kernel (fT2|) possesses an interesting dependence 
on the regularization procedure of the running coupling. This is related to the fact that this kernel is a complicated, 
non-linear function of a s . In particular, the factors of the coupling in the denominators in expression ([12]) lead to a 
non-trivial dependence on the regularization parameter of the strong coupling. It was found that by increasing the 
value of fi, see Eq. (|16p . and thus decreasing the maximal value at which the coupling freezes, there was a region of 
dipole sizes where the solution obtained from evolution with the kernel (TT21 was actually increased, contrary to what 
could be naively expected. We stress that this behavior was observed for some range of dipole sizes only. With the 
minimal dipole size prescription the amplitude was of course always decreasing with increasing value of [i as expected. 

In general, it was found that the solution with the running coupling in the form (|12p possessed a larger sensitivity 
to the way the coupling is regulated than the evolution with the minimal dipole prescription (|14[) . This sensitivity 
persists even at small dipole sizes which are far away from the scale -j. It suggests that the terms with inverse 
coupling in kernel ([T2"j) increase the sensitivity to scales which are different than the scale set by the parent dipole Xoi . 
It is also important to note that this behavior was found both in the evolution with and without impact parameter 
dependence. The solution using the minimum dipole size prescription does not exhibit such a large sensitivity. It 
would be interesting to investigate these features further in order to determine whether this is a physical behavior or 
it is an artefact of the truncation of the resummation of perturbative series which lead to this result (70j . 

Finally, the dependence of the amplitude on impact parameter for fixed value of the dipole size was analyzed. The 
diffusion property of the solution in impact parameter space is illustrated in Fig. [7] for the running coupling case in 
scenario (TT5)) . Plots in Fig 7(a) and Fig. |7(b)| differ only by the choice of horizontal scale. The solution to the BK 
equation is compared with the profile in impact parameter using the model © with the parameters from [45|. We 
observe that for the small values of b in general the BK solution is fairly close numerically to model ©. There is 
however a significant difference in the shape of the amplitudes between the BK calculation and the Glauber-Mueller 
model ([5]). This is especially manifest at large values of impact parameter where the BK solution has a more extended 
tail in b. In the BK solution there is a clear increase of the width of the distribution in impact parameter with 
increasing rapidity. 

The diffusion property in impact parameter is best illustrated in Fig. [3] where we show the average width squared, 
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Dipole Size: 1.00 (GeV"') | cos(phi): 0.0 | Delta Y: 4.0 | max Y: 8.0 



Dipole Size: 1 .00 (GeV ) | cos(phi): 0.0 | Delta Y: 4.0 | max Y: 8.0 




Impact Parameter (GeV" 1 ) 
(a) Dipole size r = 1.0 GeV -1 . Logarithmic horizontal axis. 




2 3 4 

Impact Parameter (GeV" 1 ) 



(b) Dipole size r = 1.0 GeV . Linear horizontal axis. 



FIG. 7: Dipole scattering amplitude as a function of the impact parameter for fixed dipole size and dipole orientation 6 = n/2. 
The solid lines represent the model (0 used in [4^]. The dashed lines correspond to the solution of the BK equation with the 
kernel (|15JL m — 0.35 GeV. The dashed - dotted line represents the initial conditions at Y — (xo = 0.01) also taken from 
model in 14511. 
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FIG. 8: The value of the average squared width {b 2 }, defined in Eq. (I17[) . as a function of rapidity for fixed value of the dipole 
size r. The solid line is the model J8| with parameters taken from [45| and the dashed line is obtained from solution to the BK 
equation with the kernel (|15p . 



defined as 



(b 2 ) 



J d 2 bb 2 N(r,b;Y) 
J d 2 hN{r,h;Y) 



(17) 



as a function of rapidity for fixed value of the dipole size r. We compared the value of (b 2 ) extracted from the solution 
to the BK equation with the value obtained from model ([8]) . The model ([8]) gives almost constant width, independent 
of rapidity, which is to be expected. On the contrary, in the case of the BK equation the width clearly increases with 
rapidity. For the rapidities considered here, we observe that it is almost a linear growth, with slightly faster increase 
at the highest values of rapidity ~ 6 — 8 along with mild dependence of the slope on the value of the dipole size. 



13 



B. Description of the F2 and Fl structure function data 



In our calculation of the F% and Fl structure functions we used the solution to the BK equation with impact 
parameter dependence and a gluon mass m as implemented in scenario (|15[) . The initial condition for the calculation 
was set as in Sec. IIII Al In addition, we imposed a cutoff on the dipole sizes to be equal to x 01 = ^ in the initial 
condition. As discussed earlier, this is necessary in order to be consistent with the cutoff placed in the evolution 
kernel. The data on the structure function correspond to the combined HI and ZEUS data [57] . and only the points 
below x = 0.01 are used here. 

In this calculation there are the following parameters: C, 770, Bq in the initial condition j5J), mass m in the kernel, 
strong coupling regulator \x and Aqcd- In principle all of these parameters can be varied to obtain the fit to the data. 
In practice however, due to very long time needed to find solution for a given set of parameters (of the order of a day 
on 32 cores), variation of all parameters is prohibitive. We have instead chosen to fix all of the parameters with the 
exception of fi and m which we varied in a very limited range. 

Since we use an initial condition which is cut at large dipole sizeSjthe data at values of x around 0.01 are underes- 
timated by our model. This is because the initial model flSJ) from [45[ was fitted to the data without any cuts. This 
indicates a rather large sensitivity of the F2 obtained from @ to the region of large dipole sizes even for moderate val- 
ues of Q 2 . This effect is well known and corresponds to the presence of the aligned jet configurations in the transverse 
structure function [z3j Iz3 • This can be seen by inspecting Eqs. (|2l5p as this expression receives large contributions 
from the endpoints z ~ 0, 1. As a result, at a given value of Q 2 the dipoles which contribute to F2 form a rather 
wide distribution in dipole size. We have verified that for the model given by ((5|) with a cut of the order of ^ and 
— , with m = 0.35 GeV the contribution to F2 at Q 2 = 4.5 GeV 2 from dipoles larger than the cut is about 30% and 
10% respectively. Hence, the F2 structure function contains significant non-perturbative contributions even at the 
moderate values of Q 2 . In order to compensate for this non-perturbative off-set one could of course move the cutoff 
m towards smaller masses. However, since the cut on the dipole sizes is strongly correlated with the profile in impact 
parameter it would result in the much larger width of the impact parameter profile, which would be inconsistent 
with the data on diffractive J/^S production. Therefore, we choose to work with the value of the cut which is more 
consistent with the number obtained from the J/^B diffractive slope. As a result, in order to compensate for the offset, 
a separate non-perturbative contribution is added which is important at low values of Q 2 , of the order of < 15 GeV 2 . 

We stress that this property stems from the fact that in the BK evolution equation the impact parameter and dipole 
size are strongly correlated. This has to be contrasted with the Glauber-Mueller like parametrization ([8]) where the 
dipole sizes and impact parameter are decoupled. The non-perturbative part originating from large dipole sizes was 
parametrized in the following form 



^2 S ° ft = ^—°o I rdr [ dz (|*l| 2 + |*t| 2 ) . (18) 
2na em J A. Jo 

The total structure function is then taken to be of the form 

ptot = pBK + F |oft ^ ( 19 ) 

where F 2 BK denotes the contribution obtained by using the solution to the BK equation (which is the perturbative 
part). In the formula (fT8"|) we assumed that the dominant part of the integral is where the dipole - proton amplitude 
is almost flat and therefore replaced the dipole cross section with the constant <tq ■ Note that the integral over dipole 
size in (|18[) is now cut from below by 1/m, and the integral extends into the large dipole regime. <jq is a constant 
that is used to fit the data at x = 0.01 and lowest bin in Q 2 . The F% contribution is slowly varying with Q 2 and it 
accounts well for the non-perturbative dipoles at low Q 2 . This procedure of adding separate contributions from small 
(perturbative) and large (non-perturbative) dipoles is similar to the one employed in Refs. [75L l76l|. 

One could argue that the non-perturbative contribution could be accounted for by including the contribution from 
the vector meson dominance model [77l - [80j . The vector meson dominance (VMD) contribution can be written as 



p(vmd) _ Q 2 _ v / m$<r v (W 2 ) Q 2 m 2 v a v (W 2 ) . ,„,. » , 



see for example [751 ] . Here m v denotes the vector meson mass and a v is the vector meson-proton cross section which 
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is a function of the energy. These cross sections can be taken to be 

cr p = a u = - {(t(tt + p) + (T(ir~p)) , (21) 

a* = <j(K+p) + a(K-p) - ± (a{n+p) + a^p)) , (22) 

where the parameterization for the energy dependence of (T(7r ± p), a(K ± p) was taken using the soft pomeron model 
from [81]. The j v terms relate to the leptonic width of the vector meson [751 ] in question and are defined by j 2 = 
™e m m v whgj-g the leptonic decay width is taken from 82] and £0 is a constant taken to be 0.7 [75j |. 

The computations for F 2 using the non-perturbative contribution of either (fT8|) or (|20| are presented in Figs. I9|10l 
As we see the curves with the VMD term undershoot the data whereas the curves with the soft term are systematically 
closer. The VMD contribution is almost negligible with the exception of the lowest Q 2 bin. We see that the two 
calculations, with soft term and VMD are close for values of Q 2 > 15 GeV 2 which means in this region the data are 
only described by the perturbative component. We also conclude that the VMD model is not sufficient as the only 
soft contribution for the description of the data in the present setup. VMD only contributes to the very low values of 
Q 2 (less than 4 GeV 2 ) and an additional soft component is needed which has a flatter Q 2 dependence. This result is 
consistent with the results of the previous analyses |76j . 

We observe that the slope of the calculations is too steep in x in all bins of Q 2 which implies that the LL evolution 
with the running coupling leads to a faster slope in In 1/x than the data indicate. This can be remedied by lowering 
the scale Aqcd from the value which we used, i.e. 0.25 GeV, and taking it as a fitting parameter. This was effectively 
done in the fits presented in |37J,|38(. We estimated that it would require Aqcd to be well below the pion mass, of the 
order of tens of MeV or so to fit the data. The fact that the LL evolution with running coupling in our scenario has 
a steeper slope than the data is not unexpected as one needs to take into account the next-to-leading corrections to 



BK equation. These have been computed in [72j , but a detailed analysis of the BK equation which includes them still 
needs to be performed. Preliminary analysis in the momentum space was recently performed (83j using the method 
of the saturation boundary 28]. The results from this analysis indicate that next-to- leading corrections to the BK 
equation, which are not due to the running coupling, are indeed substantial and can lead to the instabilities of the 
evolution despite the presence of saturation [84j . This strongly s ugg ests that a resummation of subleading corrections 
in In 1/x is needed in addition to the saturation corrections [85j, |86| . 

Finally, we also compared the calculation to the experimental data on the structure function Fl [56j . this is 
illustrated in Fig.[TTJ We see that the calculation is consistent with the experimental data in all bins of Q 2 . However, 
the data on Fl have very large errors. Note that, in the figures presented, the range in x in each of the bins is very 
small, and therefore the Fl structure function is very flat in each of these bins. 



V. CONCLUSIONS 



In this paper we have analyzed the nonlinear BK equation with impact parameter dependence and running coupling 
in the presence of a mass scale, which regulates dipole splitting in the infrared regime. This effective gluon mass is 
responsible for the non-perturbative effect of confinement. Using the resulting solution for the dipole scattering 
amplitude we have performed a comparison with experimental DIS data on the structure functions F2 and Fl- Let 
us summarize the main points of this investigation : 

1. The details of the evolution in rapidity strongly depend on the way the large dipoles are regularized. In particular, 
the speed of the evolution with rapidity is affected by the choice of regularization. Two different scenarios have 
been tested: modified Bcssel functions (fTU)) , and a more stringent cutoff with theta functions (fTTj) . The first 
scenario possesses a physical motivation and can be derived from the computation of the branching kernel for 
dipoles in the presence of the effective gluon mass in the propagators. 

2. The scenario (TTU1) does not entirely tame the evolution into the large dipole size region and in the presence of 
the running coupling it results in a rather fast evolution in rapidity. The resulting amplitude is also much larger 
than in scenario (fTTj) , not only in the large dipole regime but also in the small dipole region as well. We found 
that the scenario with the cutoff on all the large dipoles in the form (|11[) results in solutions which are preferred 
by the experimental data. 

3. The running coupling prescription (1121) gives, in general, much slower evolution than the prescription (|14[) with 
the minimal dipole as the scale of the running coupling. This happens despite the fact that the two kernels are 
formally equivalent in the limits when the dipole sizes are strongly ordered. We also found that the evolution 
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with kernel (fl~2]) is more dependent on the details of the regularization of the running coupling. Most likely 
this is caused by the highly nonlinear form of (TT2"1) . i.e. by the fact that it contains the inverse powers of the 
strong coupling. We found that the scenario which gives results closest to the data is the kernel with Balitsky 
prescription for the running coupling and the massive regulator taken in the form (|15p. 

4. Comparison with the data on structure function F2 shows that the slope in x of the calculation is too steep 
for the data in the case of the LL evolution with running coupling. This could be cured by changing the value 
of Aqqd, but it would require unrealistically small values of this scale in order to match the data. We stress 
that the full fit which would include the variation of all parameters, in the case of the BK equation with impact 
parameter, would be extremely demanding as far as computing resources are concerned (using similar techniques 
as presented here). The fact that the scale in the running coupling has to be adjusted quite a bit to fit the data 
is consistent with previous calculations (which were done without impact parameter dependence) and calls for 
the inclusion of the remaining corrections beyond the leading - logarithmic order. 

5. An important feature of the solution to the BK equation is the fact that the impact parameter and dipole size 
are strongly correlated. This property stems from the basic form of the dipole evolution kernel. This has to 
be contrasted with the models previously used in the literature, for example of the form l[5]). This correlation 
introduces novel effects and leads to more constraints on the calculations. For example, incorporating the 
effective gluon mass in the kernel introduces a scale in impact parameter. Strong variation of this scale is not 
possible if one requires consistency with the observed slope in diffractive data. This scale then results in the 
truncation of the large dipole size region, which in turn causes the offset in the F2 calculation, particularly 
at small values of Q 2 . Therefore one needs to include an additional (soft) component to F2 which is entirely 
non-perturbative. The BK solution also exhibits the diffusion in impact parameter, a feature that is completely 
absent in the models of the form (JSJ) . We found that in the region of large impact parameters and high rapidities 
the differences between the BK solution and the model (SJ were substantial. 
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FIG. 9: Fi proton structure function versus x in bins of Q 2 . The dipole amplitude is obtained from the BK evolution with LL 
kernel with running coupling in scenario (|15[) . fi = 0.52 GeV, m = 0.35 GeV, ao = 75.98 GeV -2 and the initial condition is 
cut at zoi = The dashed line corresponds to the calculation with the VMD term in Fi and the solid line corresponds to 
the calculation (| 181 19p . The data points are taken from the combined HI and ZEUS data sets [57) |. 
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FIG. 10: F2 proton structure function versus x in bins of Q 2 . The dipole amplitude is obtained from the BK evolution with 
LL kernel with running coupling in scenario ()15|l . /1 = 0.52 GeV, m = 0.35 GeV, ao = 75.98 GeV -2 and the initial condition 
is cut at xoi = £7. The dashed line corresponds to the calculation with the VMD term in F2 and the solid line corresponds to 
the calculation (|18I19[I . The data points are taken from the combined HI and ZEUS data sets [571 ]. 
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FIG. 11: Ft proton structure function versus x in bins of Q 2 . The dipole amplitude is obtained from the BK evolution with 
LL kernel with running coupling in scenario (|15[1 . The parameters are: \i = 0.52 GeV, m = 0.35 GeV. The initial condition is 
cut at xoi = ^r. The non-perturbative contributions are not included as they are much smaller for the longitudinal structure 
function. Note the horizontal axis range changes from bin to bin and in generally contains very small range in x. 



